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SUSY SIGNATURES IN ATLAS AT LHC 



FRANK E. PAIGE 
Brookhaven National Laboratory, Upton, NY 11973 

This talk summarizes work by the ATLAS Collaboration at the CERN Large 
Hadron Collider on the search SUSY particles and Higgs bosons and on possi- 
ble measurements of their properties. 

1 Introduction 

It has been twenty years since Richard Arnowitt, AH Chamseddine, and Pran 
Nath introduced minimal supergravity (mSUGRA) as a phenomenologically 
viable model of SUSY breaking^ This talk summarizes results from the 
ATLAS' Detector and Physics Performance TDR and more recent work by 
the ATLAS Collaboration on the search for and possible measurements of 
SUSY particles at the LHC. It also discusses measurements of Higgs bosons, 
which are a necessary part of SUSY. Much of the ATLAS work continues to 
be based on the mSUGRA model commemorated at this meeting. 

If SUSY exists at the TeV scale, then gluinos and squarks will be copi- 
ously produced at the LHC. Their production cross sections are comparable 
to the jet cross section as the same Q^; if i? parity is conserved, they have 
distinctive decays into jets, leptons, and the invisible lightest SUSY particle 
(LSP) Xi, which gives ^t- Since ATLAS (and CMS) are designed to detect 
all of these, simple cuts can separate SUSY events from the Standard Model 
(SM) background. The main problem at the LHC is not to discover SUSY but 
to make precise measurements to determine the masses and other properties 
of SUSY particles. This wiU help to understand how SUSY is broken. SUSY 
models in which R parity is violated have also been studied, but they will 
not be discussed here. 

Since the main background for SUSY is SUSY, ATLAS has emphasized 
studies of specific SUSY model points. Most of these studies start by generat- 
ing the signal and the potential SM backgrounds using a parton shower Monte 
Carlo (HerwigS, Isaje^, or Pythia^). The detector response is simulated us- 
ing parameterized resolutions and acceptances derived from GEANlEl, and an 
analysis is developed to isolate specific SUSY channels. Recently some work 
has been done using full GEANT simulation and reconstruction directly. 

2 Search for SUSY Particles at LHC 

For masses in the TeV range SUSY production at the LHC is dominated by g 
and q. Leptonic decays may or may not be large, but jets and are always 
produced, and these generally give the best reach. Consider an mSUGRA with 
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Figure 1. M^g distribution for a typical mSUGRA point and SM backgrounds after cuts. 
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Figure 2. Search limits for the mSUGRA model in various channels for lOfb"^ (left) and 
overall for various luminosities (right) 



mo = 100 GeV, mi/2 = 300 GeV, Aq = 0, tan/? = 10, sgn^ = +. Require 
> 100 GeV, at least four jets with Et > 100, 50, 50, 50 GcV, and plot as 
a measure of the hardness of the collision Afcff — + J2j ^tj ■ Then as 
Figure n shows the SUSY signal dominates for large Moff. The search limits 
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Figure 3. Search limits for the AMSB modell^l 



from this sort of analysis for mSUGRA requiring S" > 10 and S/ \fB > 5 reach 
more than ITeV for only 0.1 fb"! and 2TeV for 10 fb"!; see FigureEl 

While the AMSB modeP^is quite different, the reach in Mg, Mq is similar: 
above 2TeV for 100 fb^^. Overall reach depends mainly on a{Mg,Mq) pro- 
vided that M^o <C Mg, Mq, so one expects similar reach in most i?-conserving 
models. This should be sufficient if SUSY is related to the naturalness of the 
electroweak scale. 

3 SUSY Particle Measurements 

If R parity is conserved, all SUSY particles decay to an invisible LSP Xi, so 
there are no mass peaks. But it is possible to identify particular decays and to 
measure their kinematic endpoints, determining combinations of masses .^^^^ 
The three-body decay ~^ Xi^~^^~ gives a dilepton endpoint at Mu = 
M^o — M^o, while X2 ~^ i^i^ Xi^'^^~ gives a triangular distribution with 
an endpoint at 

M« = ^(M| - M?){M? - M^^^/M^ . 

These endpoints can be measured by requiring two isolated leptons in addition 
to multijet and f^r cuts like those described above. If lepton flavors are 
separately conserved, then contributions from two independent decays cancel 
in the combination e+e^ + /i^/i^ — e*/i^ after acceptance corrections. The 
resulting distributions after cuts. Figure 01 are very clean and allow a precise 
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Figure 5. Distributions of the smaller M (££j) (left) and the larger M(llj) for M^i > 
Mf^-^l^ (right) mSUGRA "Point 5. 'El 



measurement of the endpoint. The shape allows one to distinguish two-body 
and three-body decays. 

Long decay chains allow more endpoint measurements. The dominant 
source of X2 mSUGRA "Point 5'^ and similar points is qi^ x"? ^ 
i^P^q y^l'^l~ q. Assume the two hardest jets in the event are those from 
the squarks and for each calculate M(e£j), M<{£j), and M>{£j). Then the 
smaller of each of these should be less than the endpoint M«g, M^^\ ^eq'^ 
for squark decay, while the larger M{Uj) should be greater than the threshold 
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Figure 6. left: Scatter plot of A/; vs. Af-o for two models consistent with measurements 
described in the text. Right: Projection of M-o 

Tegq requiring Mu > Mf^'^^ These endpoints are smeared by jet recon- 
struction, hadronic resolution, and mis-assignment of the jets that come from 
squark decays. Nevertheless, the distributions show clear structure at about 
the right positions. 

After accumulating high statistics and careful study, it should be possible 
to measure the endpoints to the expected hadronic scale accuracy, ^ 1%. The 
Uq threshold is more sensitive to hard gluon radiation, so it is assigned a larger 
error, ^ 2%. Some distributions of the resulting masses derived assuming 
these errors are shown in FigureElfor two models, mSUGRA Point 5 (S5) and 
an Optimized String Model (01) with similar similar masses. Relations among 
the masses are determined to ~ 1% and are clearly sufficient to distinguish 
these models. The LSP mass is determined to ~ 10% by this analysis; since it 
is determined only by its effect on the kinematics of the decay, the fractional 
error on M^o clearly diverges as M^o/M~ — > 0. 

4 h ^ bb Signatures 

If X2 ~^ Xi^ is allowed, it may dominate over X2 ~^ Xi^^- This signal can 
be reconstructed using two b jets measured in the calorimeter and tagged as 
6's with the vertex detector. A typical signal using the expected 6-tagging 
efficiency and light-quark rejection and the reach for such signals are shown 
in Figure Ul Such a Higgs signal in SUSY events might well be observed with 
less luminosity than h — > 77 or h ZZ* and so be the discovery channel for 
the light Higgs. 

If a signal for h ^ bb is observed in SUSY events, the h can be combined 
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Figure 8. Mf^j distribution for events with h — > bb\ 



with the two hardest jets in the event to measure the q —>■ x\hq endpoint in a 
way similar to the measurement of the U.q endpoint. The resulting distribution 
is shown in Figure |H1 the endpoint is consistent with what is expected. While 
the errors are worse than for the U,q endpoint, the measurement is still useful. 
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5 Heavy Gaugino Signatures 

In mSUGRA and other typical SUSY models, the light charginos and neu- 
tralinos are mainly gaugino and so dominate the cascade decays, so that 

B{qL ^ X°2l) - 1/3, BiQL ^ xfq') - 2/3, B{qn ^ x??) - 1 ■ 

But even in the simplest mSUGRA model, Xi ^-nd X2 have a significant 
admixture of gaugino and so contribute in light-quark decays of squarks and 
gluinos. 

Four xl/xf decay chains can give OS, SF dileptons: ql — > XaQ 
ii^^q xlt+^-q [Dl]; qLxlq ^ i^i^q ^ x\i+ q [D2]; uxlq ^ i^i^q ^ 
xl^+^-q [D3]; and qL X2I' ^ ''e^^l' ^ Xi^'^q' [D4]. In principle these 
four decay chains give four distinct £^£^ endpoints, but it seems impossible 
to resolve these even with 100 fb^^ of integrated luminosity. Nevertheless, 
there are > 10'^ events from heavy gauginos over substantial range of 

mSUGRA parameters; Xi decays dominate for low mo, while X2 dominates 
for the region mo '■-^ rni/2- 

Event samples were generated and simulated for each of the ten points 
indicated in Figure El Events were required to have an £~^£~ dilepton pair, 
Mu > 100 GeV, > 100 GeV, > 4 jets, and M^g > 600 GeV. To suppress 
SM backgrounds a cut Mt2 > 80 GeV was also made, wher^^ 

M^2 = min [ma.x{mTipT4i^M^'i^TiPT,i2^h)}] 

is the minimum transverse mass obtained by partitioning the observed fir 
between two massless particles. Note that Mt2 < Mw for t and W back- 
grounds. 
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Results for mo, toi/2 = 100, 250 GeV and 150, 250 GeV are shown in Fig- 
ure ^1 Evidently the signal is observable over the SUSY and SM background 
in both cases. The estimated statistical error on the cndpoint is about ±4 GeV 
in both cases. The 5a reach for such signals in mSUGRA is indicated by the 
dark curve in Figure ITTl Heavy gaugino signals are rather model dependent, 
so the ability to study them is important for understanding the SUSY model. 
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Figure 12. Rec onst ructed Mti, distribution before sideband subtraction (left) and after sub- 
traction (right) EI 



6 Third-Generation Squark Signatures 

The properties of the third-generation squarks 61,2 and ti_2 are important for 
understanding the SUSY model, but their signatures are typically complex. 
The main production mechanism is g production and decay. Consider for 
mSUGRA with mo = 100 GeV, mi/2 = 300 GeV, Aq = -300 GeV, tan/3 = 
10, sgn At = + the processes 

Then the M{tb) endpoint can be used to measure a combination of masses of 
the squark mass es. 

The analysiJ2l requires as usual multiple hard jets and large If^T plus two 
jets tagged as 5's and two other jets j not tagged as 5's and consistent with 
tb jjbb. The resulting Mj^ distribution is still dominated by combinatorial 
background. The next step is to select sidebands around Aljj = Mw, rescale 
the jet momenta to Mw, and subtract to determine tb signal. The M{tb) 
mass distributions for one point before and after subtraction are shown in 
Figure IT^ The fitted endpoint for this case is 443.2 ± 7.4 GeV compared to 
expected 459 GeV. A similar agreement between reconstructed and expected 
endpoints was found for all twelve points studied. Heavy squark signatures 
are clearly difficult, but it appears possible to use a sideband analysis such as 
this to study them with the ATLAS detector. 

7 T Signatures 

The mSUGRA model assumes e-/2 universality, and this is certainly suggested 
by the stringent limits on — )■ 67. Even in the simplest mSUGRA model, 
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Figure 13. t identification in ATLAS 



however, the f behave differently than e and ji because of Yukawa contri- 
butions to the RGE's, gaugino-Higgino mixing, and f]^-fj^ mixing, which is 
cx rrir- Hence r's provide unique information and might even be dominant in 
SUSY decays. 

The ATLAS (and CMS) vertex detectors cannot cleanly identify r Ivv, 
so it is necessary to rely on hadronic r decays. The background for such 
decays is much larger than that for electrons and muons. The r efficiency 
vs. jet rejection shown in Figure El should be compared with the > 10^ 
rejection for 90% efficiency expected for electrons and muons. Furthermore, 
all r decays contain missing neutrinos. For A tt one can project fix on 
the measured r directions to reconstruct the rr mass, but this is not possible 
for SUSY because of the dominant f^x from the Xi's- 

Decays into r's are generally enhanced for tan/3 ^ 1. A mSUGRA 
model with tuq = mi/2 = 200 GeV, Aq = 0, tan/3 = 45 gives X2 ~^ t'it ^nd 
xf — > tiVt with branching ratios close to unity. For events from this point, a 
simple model for the detector response turns a sharp edge at Mrr = 59.64 GeV 
into the distribution shown in Figure 1141 The visible momentum or mass 
depends both on the momentum and on the polarization of the r. Measuring 
the r polarization requires separating different r decay modes; the visible 
energy depends strongly on polarization for r — > tti/ but weakly for t aiv. 
Such a separation of decay appears to be possible, albeit difficult: for example 
T —^ TTv has a single track with p — E and low electromagnetic energy. Recent 
work based on full GEANT simulation has given an encouraging indication 
that the rr endpoint can be inferred from the visible tt mass. 
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Figure 14. Reconstructed visible tt mass in a mSUGRA event sample 



8 GMSB Signatures 

While the mSUGRA model remains after 20 years perhaps the most attrac- 
tive paradigm for SUSY breaking, it may not be correct. In the GMSB 
model SUSY breaking is communicated via gauge interactions at a scale 
much less than the Planck scale, so the gravitino G is very light. GMSB 
phenomenology'^ depends on the nature and lifetime of NLSP (x? or I) to 
decay into the G. In general the GMSB model produces longer decay chains 
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with more precisely measured decay products, so reconstructing masses is 
considerably easier than in mSUGRA. 

The GMSB model can give a number of special signatures related to 
NLSP decay. If the NLSP is a x^i its lifetime for Xi ~^ can range from 
very short to very long. Short lifetimes can be detected by using the Dalitz 
decays x? G£'^ell~ with branching ratios of a few percent. Long lifetimes 
can be detected by looking for (rare) non-pointing photons in SUSY events. 
The ATLAS electromagnetic calorimeter has both good angular resolution in 
the polar angle, A6 ~ , ^° , and good timing resolution, At « 100 ps. 

f t> 1 y^E/lGoV b b T f 

Both can be used to detect non-prompt photons from long-lived particles like 
Xi produced with /3 < 1. Such signals give a sensitivity up to ct ~ 100 km, 
much greater than what is expected in the GMSB model. 

For other choices of the parameters the GMSB model might give long- 
lived sleptons, which look like muons with /3 < 1 in a detector. The ATLAS 
muons chambers give a timc-of-flight resolution in the 1 ns range over a dis- 
tance of about 10 m, making it possible to reconstruct both the momentum 
and the mass of the slepton. The slcpton lifetime can be determined by 
comparing the rates for events with one and two reconstructed sleptons as 
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shown in Figure 1161 The statistical error is small; the dominant systematic 
error is difficult to estimate without real data. Another approach would be 
to look for sleptons decaying into non-pointing tracks in the central detector. 
This should be more sensitive for long lifetimes, but estimating the sensitivity 
requires studying the pattern recognition for such non-pointing tracks. 
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Figure 17. e+e 



mass distribution from full GEANT simulation after 



corrections for average acceptance and electron energy scalell2l 



9 Full Simulation of SUSY Events 

Most studies of SUSY signatures in ATLAS have been based on fast simulation 
such as ATLFAST. While this should represent the ultimate performance of 
the detector, it does not necessarily represent the effort needed to achieve 
that performance. Therefore, a sample of 100k SUSY events has recently 
been simulated with full GEANT for an mSUGRA point with 

Too = 100 GeV, mi/2 = 300 GeV, = -300 GeV, tan/3 = 6, sgn/i = + 

The simulation of each event takes about 10'^ s, compared with about Is for 
event generation and fast simulation. Thus, such a study represents a large 
effort. 

Most of the effort so far has been devoted to debugging the reconstruction 
software, so the results are not yet useful for assessing the performance of the 
ATLAS detector for SUSY. However, a few physics plots have been produced 
using cuts based on previous fast simulation studies like those described above. 
As an example. Figure [T7I shows the e^e~ + n'^ n~ — e^/i^ mass distribution 
from full simulation of SUSY events after corrections for the average e and 
acceptance and for the e energy scale. It is encouraging that the distribution 
is quite similar to that obtained from fast simulation. 

10 Higgs Signatures 

SUSY requires Higgs bosons. In the Minimal Supcrsymmctric Standard Model 
(MSSM) there are two Higgs doubles and hence after electroweak symmetry 
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Figure 18. Significance of SM Higgs modes in ATLAS for 100fb-l|2l 



breaking five Higgs bosons, h, H, A, and H^. Tlie light, CP-even, h satisfies 
Mh < Mz at tree level Alh ^ 130 GeV after loop corrections. In many 
although not all SUSY models the h is very similar to a SM Higgs of the same 
mass. 

The search for SM-like Higgs bosons has been a principle design goal 
of both ATLAS and CMS, and a large amount of effort has been devoted 
to studies of how to search for such particles. The global summary of these 
studies is shown in Figure El for each SM Higgs mass there is at least one 
channel giving a significance of more than 5(7 for an integrated luminosity of 
100 fb^^, and the combined significance of all channels is greater than about 
lOcr. 

Recently, more effort has been devoted to studies of how to measure the 
properties of Higgs bosons once they are discovered. The key for doing this 
is to observe the Higgs boson in more than one production and/or decay 
channel.- While gg —+ h is the dominant production process at the LHC, 
WW ^ h is also significant and plays a crucial role in the analysis. These 
events can be identified by requiring hard forward jets resulting from the 
radiation of the W^s from incoming quarks, q — > Wq', and no additional 



14 




central jets. 

A typical result from such a study is shown in Figure 1191 Events were 
selected using a combination of e/^, ee + fifi, and £h modes requiring a double 
forward jet tag and a central jet veto; Mtt was reconstructed by projecting 

on the measured r directions. The accepted cross section is about 1.0 fb 
on a total SM background of 0.5 fb. Thus this channel can be used to measure 
the product of Th^ww and T^.tt- Combining a number of such measurements 
can give a good determination of the properties of the h, although a linear 
collider with sufficient energy and luminosity could do better. 

It may also be possible to reconstruct heavy Higgs bosons, especially for 
tan/3 ^ 1. An example of the reconstruction of gb — + H^t with 
T^v is shown in Figure TK\ Since the signature is a narrow hadronic r plus 
reconstructed events were selected requiring an identified hadronic r jet 
plus two non-6 jets and a h jet consistent with t kinematics. This analysis 
relies on the fact that in tuv tt^vv the vr^ is hard and so is well 

separated from SM backgrounds. 



11 Outlook 

If SUSY exists at the TeV mass scale, ATLAS should find signatures for it 
quite easily at the LHC. If R parity is conserved, no mass peaks for SUSY 
particles can be reconstructed, but several techniques have been developed 
to measure combinations of SUSY masses using kinematic distributions of 
observable decay products. 

While the details of SUSY analyses at the LHC certainly depend on the 
details of the SUSY model, it is possible to sketch a general outline of how 
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ATLAS could proceed first to search for and then to study SUSY with im- 
parity conservation: 

1. Search for an excess of niultijct + events over the SM expectation, 
and observe the Mcs at which this emerges from the background. 

2. If such an excess is found, select a SUSY-dominated sample using simple 
kinematic cuts. 

3. Look in this sample for special features such as prompt 7's or long-lived 
£; either of these may occur in GMSB. 

4. Look in the SUSY-dominated sample for £^£^, b jets, hadronic 
t's, etc. 

5. Try simple endpoint-type analyses. 

Carrying out such an initial study seems quite feasible. Its results would of 
course guide further more detailed and more model dependent analyses. 

I thank my many ATLAS collaborators who have contributed to the 
work presented here. This work was supported in part by the United States 
Department of Energy under Contract DE-AC02-98CH10886. 
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